Journal of Biological Physics manuscript No. 

(will be inserted by the editor) 



Damian H. Zanette • Sebastian Risau Gusman 



> 

oo 
o 



o 



Infection spreading in a population with evolving 



^ : contacts 
o 

> 
O 



Received: date / Accepted: date 



Abstract We study the spreading of an infection within an SIS epidemiological model on a network. 
P I , Susceptible agents are given the opportunity of breaking their links with infected agents. Broken links 

are either permanently removed or reconnected with the rest of the population. Thus, the network 
coevolves with the population as the infection progresses. We show that a moderate reconnection 
. frequency is enough to completely suppress the infection. A partial, rather weak isolation of infected 

agents suffices to eliminate the endemic state. 
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1 Introduction 



Outbursts of epidemics in human populations trigger individual and collective reactions that can 
substantially alter the social structure. As a consequence of risk perception, non-infected individuals 
may start avoiding contact with their infected equals, even when their previous relationship was fluent. 
The whole society could collectively decide to isolate its infected members until danger is overcome. 
More altruistic non-infected individuals may be tolerant of the contact with infected individuals but, in 
turn, the latter may discontinue the relationship to impede contagion. The escape from crowded cities 
during the Black Death in the late Middle Ages, documented in Giovanni Boccaccio's Decameron, and 
the closing of schools, churches, and theaters during influenza epidemics in the early twentieth century, 
constitute dramatic historical instances of such behaviours [HE]- Quarantine protocols, and preventive 
isolation during leprosy or tuberculosis treatment, are present-day examples [3]. In any case, these 
changes in the pattern of social contacts help to limit and control the incidence of the infection. 

In this paper, we explore the effects of an evolving pattern of contacts on the dynamics of infection 
spreading, in the framework of a simple epidemiological model. We consider a population of agents 
whose pattern of contacts is represented by a network. If a link of the network joins two agents, 
contagion is possible when one of them is infected and his neighbour is not. To account for the social 
processes addressed in the previous paragraph, we admit that the contact network is not a static 
structure, but evolves in response to the epidemiological state of the population. 

Agent-based models whose interaction patterns are represented by networks have received increasing 
attention during the last years, in the analysis of emergent collective behaviour in complex systems 
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[4]. Frequently, the evolution of the interaction network and the dynamics of individual agents occur 
over different time scales. In learning processes, for instance, connections change adaptively over scales 
that are large as compared with the internal dynamics of agents 5 . At the opposite limit, in models of 
network growth, the pattern evolves in the absence of any dynamics related to the agents pQ. When, 
on the other hand, the dynamical time scales of a population of agents and its interaction network are 
comparable, we can speak about their coevolution [6,7,8, 9]J10| . In this context, the model considered in 
the present study can be regarded as an illustration of the coevolution of agents and networks, inspired 
in the dynamics of infection spreading. 

Our model is based on an SIS epidemiological process where, at a given time, each agent can be 
susceptible (S) or infected (I). In the standard SIS process, each I-agent spontaneously recovers and 
becomes susceptible at a fixed rate, say, with probability 7 per unit time. An S-agent, in turn, becomes 
infected by contagion from his infected neighbours. If the contagion probability per unit time and per 
infected neighbour is p, an S-agent with ki infected neighbours becomes itself infected with probability 

pidt = 1 - (1 - p dt) kl -> kjp dt (1) 

during the interval dt. Within a mean-field description, if the average number of (both S and I) 
neighbours per agent is k and the fraction of I- agents is n\, we have fci = ktij. The mean- field evolution 
equation for n\ thus reads 

ii\ = -7TI1 + Anin s , (2) 
where rig = 1 — n\ is the fraction of S-agents, and A = kp. In this description, for asymptotically long 
times, n\ vanishes if A < 7. Therefore, the infection is suppressed as time elapses. If, on the other hand, 
A > 7, the fraction of infected agents approaches a finite value njf = 1 — 7/A > 0, and the infection is 
endemic. The transition between these two regimes occurs through a transcritical bifurcation at A = 7. 

In the following, we complement the standard SIS model with the possibility that the network 
of contacts changes in response to the infection spreading. Specifically, links between susceptible and 
infected agents can be broken, and either removed or reconnected to other agents. As expected, we 
find that this mechanism decreases the infection level, and can eventually suppress the endemic state. 
With respect to the standard model, however, infection suppression for high rates of contact change 
occurs through a tangent bifurcation, which in turn gives rise to a bistability regime. In this regime, the 
infection persists or dies out depending on the initial fraction of infected agents. More unexpectedly, 
infection suppression does not require a drastic overall change in the network structure, but is reached 
with a moderate unbalance between the mean number of neighbours of susceptible and infected agents. 
As we discuss in the final section, these features are robust under several variations of the dynamical 
rules. 

2 Evolution of the network of contacts: Link removal 

We address first the case where before the interaction between an S-agent and an I-agent joined 
by a network link effectively takes place, so that contagion becomes possible, they may decide to 
definitively delete that link, thus avoiding any further contact. In this situation, as far as the number 
of I-agents does not vanish, the network of contacts keeps loosing its links. According to our discussion 
of the standard SIS model ([2]), however, we expect that when the number of neighbours per agent has 
decreased sufficiently, the infection dies out. Once no I-agent remains in the population, removal of 
links stops and the systems reaches a static, fully healthy state. To quantify the mechanism of link 
removal, we assume that each link between an S-agent and and I-agent is deleted with probability q 
per time unit. 

Our formulation of the dynamics of the present model is analogous to the mean-field approach 
used to derive Eq. ((2]). We assume that the population consists of N agents, and call N\ and Ng the 
number of I and S-agents, respectively, so that N = Nj + N$. Additionally, we must introduce new 
variables to describe the structural state of the network and its relation with the epidemiological state 
of the population. Therefore, we consider the quantities Mn, Mig, and M$s, respectively, the number 
of network links joining two infected agents (II), an infected agent and a susceptible agent (IS), and 
two susceptible agents (SS). The total number of links is M — Mn + Mis + Afss- Due to removal of 
IS-links, the average change in Mis P er time unit is 



Mis = —qMjs- (link removal) 



(3) 
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The processes of recovery and infection are the same as discussed for the standard SIS model (J2J> in 
the Introduction. Recovery of an I-agent occurs with probability 7 per unit time and, for each S-agent, 
the infection probability per unit time and per infected neighbour is p. In a recovery event, when an 

I- agent becomes susceptible, there is not only a decay in the number of I-agents, but also a change in 
M11, Mis, an d Mss- In fact, the links joining the recovered agent with I and S-agents pass, respectively, 
from the II-type to the IS-type, and from the IS-type to the SS-type. The number of links of each type 
associated to a given agent is calculated using mean-field-like averages. For instance, the number of 

II- links associated to an I-agent is estimated as 2Mn/Ni, Similarly, the number of IS-links associated 
to an I-agent is Mis/Ni. Using this kind of arguments, we obtain, for each variable, the average change 
per time unit due to recovery: 



Ni = -N s = -jN h 
Mu = -27M11, 
Mis = 27M„- 7 Mi S , 
M ss = j Mis. 



(recovery) (4) 



Note that Mu + Mis + Mss = 0, because recovery events do not change the number of network links. 

To calculate the contribution of infection events in the mean-field approximation, we must evaluate 
the average number of infected neighbours of a susceptible agent. For a randomly chosen S-agent, this 
number is given by the ratio Mis/Ns- However, it should be taken into account that, to become infected 
by contagion, a susceptible agent must have at least one infected neighbour. This would restrict the 
calculation of the number of infected neighbours to those S-agents with at least one IS-link. For the 
sake of simplicity, we shall still estimate the average number of infected neighbours per S-agent as the 
above ratio, with the proviso that the approximation is valid when the overall number of I-agents is 
not too small, so that contagion is in principle possible for all S-agents. The change per time unit due 
to infection for each variable turns out to be 

Ni = -N s =pMi S , 

Mii=pMyN s , (infection) (5) 

Mis = P(2M SS - M m )Mis/N s , ^tectum) ^ 

M ss = -2pM ss Mi S /iVs. 

Again, M u + M IS + M ss = 0. 

To obtain differential equations of the type of Eq. @, it is convenient to define the fractions 
rii = Ni/N and my = My/M, with = {I, S}. In calculating the variation of my per time unit, 

we must take into account that also the total number of links M varies with time: 

Ma Ma • . . 

mu = — - 1 4-M. (6 

3 MM 2 y ' 

Since the total number of links changes by the removal of IS-links only, we have M = —qMis- On the 
other hand, the total number of agents TV remains constant. Also, by definition, we have ni + ns = I 
and mu + mis + mss = 1, so that we can limit ourselves to study the evolution of ni, mu, and mis- 
The evolution equations resulting from the above considerations are 

n[ = —ni + Amis, 

mjj = — 2mn + AmT + qmumis, (7) 
m' ls = 2m n - (1 + g)mi S + Am7 1 mi S (2m S s - mis) + <7mf s , 

where primes indicate differentiation with respect to the rescaled time t' = 7 t. We have also defined 

q = q/ 1 , A = fcp/27, (8) 

where k = 2M/N is the overall mean number of neighbours per agent. Note that, since the number of 
links M varies with time, the coefficient A is itself time-dependent. Its evolution is determined by the 
variation of M, and obeys 

A' = -gmisA. (9) 
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The numerical solution of Eqs. ([7|) and ([9]) confirms the expectation that the infection dies out as a 
consequence of the sustained removal of links. The fraction of infected agents asymptotically vanishes 
with time and, accordingly, m\i and mis also tend to zero. Meanwhile, the total number of links 
approaches a constant Mr. The number of remaining links Mr depends both on the initial number of 
links, Mo, and on the initial fraction of infected agents. In fact, for a given rate of link removal, the 
suppression of a higher initial infection level is expected to take longer times and require more removed 
links. 

Figure Q] illustrates the dependence of the remaining fraction of links, Mr/Mo, on the normalized 
rate of link removal q and on the infectivity. The infectivity is here characterized by the initial value 
of A, Ao = fcop/27, with fco = 2Mq/N. The results correspond to an initial condition with a fully 
infected population, n,j(0) = 1, and a fully connected network, Mo = N(N — l)/2, so that ran = 1 and 
mis = ™ss = 0. As expected, the fraction of remaining links decreases both with Ao and q. Higher 
infectivities require that a larger fraction of contacts is deleted before the infection dies out, and a 
larger removal probability contributes in the same direction. 

From the viewpoint of the dynamics, the case considered so far -where links keep being removed 
as long as a fraction of the population remains infected- is not especially interesting. In particular, 
the transition between infection suppression and persistence observed for q = when the infectivity 
grows, disappears as soon as the removal probability becomes positive. In a real population, moreover, 
it is expected that discontinued contacts are replaced, at least to some extent, by new social links, in 
such a way that the structure of society is not too much deteriorated. In the next section, we study a 
model where links between susceptible and infected agents are not permanently removed, but rather 
reconnected to other members of the population. Under these conditions, both infection suppression 
and the endemic state are possible, and can even coexist for a given set of parameters. The dynamics 
is accordingly richer, and new critical phenomena separating both regimes appear. 
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3 Reconnection of links 

We consider now that, before each susceptible-infected interaction -possibly leading to contagion- 
takes place, the susceptible agent is given the opportunity of breaking the contact with his infected 
neighbour and to reconnect the corresponding link to another agent, randomly chosen from the rest of 
the population. Reconnection of each IS-link occurs with probability r per time unit. If the susceptible 
agent is reconnected to another susceptible agent, the link changes from the IS-type to the SS-type. 
Otherwise, no change occurs. 

A similar model was considered recently |12j where, however, reconnection of S-agents always 
occurs towards other S-agents. In this model, upon reconnection, IS-links always change to SS-links. 
This variant is implicitly admitting that agents have information on the (S or I) state of their equals 
before making contact, which seems to be a rather artificial assumption. On the other hand, in our 
model reconnection is done at random, which overcomes such assumption but, at the same time, limits 
the efficiency of I-agent isolation. 

Only the variables Mis an d Mss change due to reconnection events. Since the probability of choosing 
an S-agent at random is ns we have, per unit time, 

Mis = — Mss = — rnsMig. (link reconnection) (10) 

Changes due to recovery and infection are the same as in Eqs. ((4]) and ([5]), respectively. 

Obtaining the evolution equations for the fractions n\, mu, and mis, is now simpler than in the 
case of link deletion, because the total number of links M remains constant. We find 

n\ = — m + Amis, 

m'u = — 2m\i + Ang 1 mj S , (11) 
m' ls = 2m u - (1 + fn s )mi S + Arig 1 mi S (2m S s - mis)- 

Here, again, primes indicate differentiation with respect to the rescaled time t' = ~ft. Moreover, 

r = r/7, A = kp/2j, (12) 

with k = 2M/N the average number of neighbours per agent. Since, now, M does not vary with time, 
both k and A are constants. Note that the normalized reconnection probability f and infectivity A are 
the only parameters of Eqs. (fTTj) . The normalized infectivity incorporates the only network-specific 
feature, namely, the mean connectivity k. 



3.1 Infection level at equilibrium 



We focus the attention on the equilibrium solutions of Eqs. (jlip . which are the candidates to represent 
the infection level and network structure at asymptotically long times. First, we consider the stationary 
values of the fraction of infected agents. The analysis is restricted to the case of A > 1/2 which, in the 
absence of reconnection events (r = 0), corresponds to an endemic infection, n\ > for t — > oo (cf. the 
discussion of the standard SIS model in the Introduction). 



At the fixed points of Eqs. (fTTj) . the equilibrium fractions of links mj 
equilibrium fraction of infected agents nj" as 



and 



'is 



are related to the 



,*2 



2A(1 



ns 



(13) 



In turn, rij satisfies 







2A - 1 - f + (3f - 2AK - 3fnf + fn\ 



(14) 



This polynomial equation has four solutions. One of them tends to infinity for f — > 0, and remains real 
and larger than one for any positive f. Since meaningful solutions to our problem must verify rij < 1, 
we disregard this solution from now on. 
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Fig. 2 Bifurcation diagram for the equilibrium fraction of infected agents n\ as a function of the normal- 
ized reconnection probability f , for three normalized infectivities A. Although only positive values of nj are 
meaningful, an interval in the negative domain is also shown for completeness. Full and dotted lines represent, 

respectively, stable and unstable branches. For clarity, the solution nj ' = is plotted in the vicinity of the 
transcritical bifurcation only. 



The trivial solution nj ' = exists for any value of the normalized infectivity A and of the normalized 
reconnection probability f. For a given A, its stability depends on f. As discussed in more detail below, 
= is unstable for small r and becomes stable as r grows. The other two solutions read 




01 -r- R ■ 01 

cos — =F v3sm — 



(15) 



with 



/32A 3 , , 

a = arctan 1/ 1 (16) 

V 27f x ' 



(0 < a < it/2). These two solutions are real for 32 A 3 > 27 f. Otherwise, they are complex conjugate 
numbers. The solution rij , with the minus sign in the right-hand side of Eq. (JTSJ), approaches 1 — (2A) _1 
for f — > 0. Thus, it represents the expected fraction of infected agents in the absence of reconnection. 
When it is real, it satisfies n.j < 1, and it is stable as long as it remains positive. Consequently, along 
with the trivial solution, rij is another meaningful equilibrium solution to our problem. Finally, nf ^ 
is negative and stable for small f. Depending on A, it can become positive as f grows but, at the same 
time, it becomes unstable. Therefore, it does not represent a meaningful solution. 

Figure [5] summarizes, in a bifurcation diagram, the behaviour of n^, and as functions of 
the normalized reconnection probability f, for three representative values of the normalized infectivity 
A. In the three cases, we have A > 1/2, so that -as discussed above- a non-trivial meaningful solution 
does exist. Full and dotted lines represent, respectively, stable and unstable branches. For small infec- 
tivity (A = 0.6), the stable solution nj 1 ^ crosses n|°^ and becomes negative and unstable, while nj ^ 
becomes stable. This transcritical bifurcation takes place at f = 2A — 1. As f grows further, rij and 
the negative stable solution nj 2 ' approach each other, and collide when nj 1 ^ = — 1 — 3/4A. Beyond 
this tangent bifurcation, which takes place at f — 32A 3 /27, the two solutions are complex numbers. 
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the tangent (Tg) bifurcation lines, which collapse into a pitchfork bifurcation (P) at (1/2, 3/4). The dotted line 
is the continuation of the tangent bifurcation line in the zone where is negative. 



Phase diagram in the (f, A)-plane, showing the regions of infection suppression (ni — > 0) and persistence 
and the intermediate bistability zone. Their boundaries are given by the transcritical (TC) and 



The situation is different for larger infection probabilities, as illustrated by Fig. [5] for A = 1. Now, 
for r = 0, 7ij is large and, as f grows, it is the stable negative solution nj 2 -* which first reaches n[ ^ . 
At the transcritical bifurcation at r = 2 A — 1, n'p becomes positive and unstable, and nf^ becomes 

stable. The tangent bifurcation where and rhp collide and become complex, at f = 32A 3 /27, takes 
now place when these two solutions are positive. As a consequence, there is an interval of normalized 
reconnection probabilities, between the two bifurcations, where the system is bistable: both rij and 

rij are stable meaningful solutions to the problem. The asymptotic state is selected by the initial 
condition for m. 

The regimes of small and large infect ivity are separated by the critical value A = 3/4 = 0.75, 
also shown in Fig. [2j At this critical point, the transcritical bifurcation and the tangent bifurcation 
collapse into a pitchfork bifurcation at f = 1/2. Here, the three equilibria collide simultaneously, and 
n[ becomes stable, while the other two solutions become complex. 

A phase diagram of our system over the parameter plane (f , A) is shown in Fig. [3J The zones 
of endemic infection, where the fraction of infected agents at asymptotically long times is positive 

{n\ — > nj 1 '), and of infection suppression (m — > 0) are separated, for large A and f, by the bistability 
region, where the two asymptotic behaviours can be obtained, depending on the initial condition. The 

three zones are limited by the lines of the transcritical bifurcation [A = (1 + f)/2, TC], where nj "* = 

changes its stability, and of the tangent bifurcation [A = (27f/32) 1 / 3 , Tg] where and rij collide 
and become complex. These two lines are tangent to each other at the "triple point" (1/2, 3/4), where 
the bistability region disappears, and the system undergoes a pitchfork bifurcation [P]. For smaller 
A and f bistability is no more possible, and the zones of infection persistence and suppression are 
separated by the transcritical line. The tangent bifurcation takes now place at negative values of 

(2) 

and nj (dotted line). 

Let us summarize our results on the persistence or suppression of the infection in terms of the non- 
normalized parameters. First, for small infcctivity, p < 7/fc, the infection is always suppressed. In this 
situation, the infectivity is just too small to sustain a finite infected population. For larger infectivities, 
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on the other hand, the infection can become established, depending on the reconnection probability r. 
In the range 7/fc < p < 3-f/2k, the infection is endemic if reconnections are infrequent, r < kp — 7. 
Otherwise, for sufficiently frequent reconnections, the infection dies out. The transition between the 
two situations is continuous in the fraction of infected agents, and occurs through a transcritical 
bifurcation. For even larger infection probabilities, p > 37/2fc, the regimes of persistence (low r) 
and suppression (large r) are separated by a bistability zone, where the infection persists or dies out 
depending on the initial fraction of infected agents. The bistability zone is limited by the transcritical 
bifurcation quoted above and a tangent bifurcation at a reconnection probability r = 4k 3 p 3 / 27 "f 2 . The 
discontinuous nature of the tangent bifurcation implies that the endemic state present in the bistability 
zone disappears abruptly at the boundary, with a finite jump in the asymptotic fraction of infected 
agents, from n\ = 1 — y/kp/ir > to zero. 



3.2 Number of neighbours of infected and susceptible agents 

The variables m\\ and mis characterize how the structure of the network is related to the state of 
the agents. Reconnection events favor the growth of the number of SS-links at the expense of IS-links. 
Thus, for r > 0, S-agents should asymptotically posses relatively large numbers of neighbours. The 
equilibrium values mjj and m| s as functions of the equilibrium fraction rij of I-agents are given by 
Eqs. (|13p . These equations show, as expected, that the fraction of links connecting I-agents with any 
other agent is proportional to the fraction of I-agents itself. 

In order to introduce quantities that define the connectivity of I-agents and S-agents independently 
of their respective fractions, we consider the average number of neighbours per agent of each type. For 
I-agents, for instance, the average numbers of infected and susceptible neighbours are 2Mu/N\ and 
Mis/Nj, respectively. The average connectivity of I-agents, fci, is the sum of these two quantities or, 
equivalently, 

h 1 

-r = — . 17 

k 2A(l-nf) 

which gives the ratio between k\ and the overall average connectivity per agent, k = 2M/N. With 
analogous arguments for S-agents, their average connectivity reads 

k s 1 /J? 



k 1 



2\{1 -n*) 2 



(18) 



Due to the conservation of the total number of links, k\ and fcg are univocally related. This relation 
can be obtained from Eqs. (fl7|) and (| 18[) by eliminating nj, which yields 

fc s = fci[l + 2A(l-fci/fc)]. (19) 

In order to describe the correlation between the structure and the state of the population it is however 
useful to analyze both k\ and ks as functions of the relevant parameters. Figure 0] illustrates the 
behaviour of k\ and fcg , as described in the following, for the infection probabilities A already considered 
in Fig. H 

For rij = nj 1 ^, which stands for the stable equilibrium solution for low reconnection probabilities, 
both ki and fcs approach k as f — > 0. As expected, in the absence of reconnection events, there is no 
difference in the number of neighbours of infected and susceptible agents. As f grows from zero, we 
have k\ < k < fcg. We thus verify that reconnection tends to increase the connectivity of S-agents at 
the expense of I-agents. 

The other solution relevant to the process, nj" = nj " 1 = 0, corresponds to a purely susceptible 
population. Accordingly, we find ks = k. Note also that Eq. (TT7|) predicts fci = k/2X, but this value is 
never realized due to the total absence of I-agents in this state. 

For A < 3/4, the fraction of I-agents decreases monotonically with f and vanishes continuously at 
the transcritical bifurcation -or, for A = 3/4, at the pitchfork bifurcation. The connectivity of S-agents 
is fcs = k both at f = and at the bifurcation. For intermediate values of the reconnection probability 
ks is larger than k and attains a maximum. This maximum, which at first glance may seem surprising, 
can be easily explained. In fact, to sustain a value of fcs larger than the overall average fc, it is necessary 
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Fig. 4 Connectivity of infected and susceptible agents, k\ and ks, relative to the overall connectivity k — 
2M/N, for three values of the normalized infectivity A, as functions of the normalized reconnection probability 
r . Only the values corresponding to meaningful stable solutions for the fraction of infected agents at plotted. 
The connectivity of infected agents is not plotted beyond the threshold of infection suppression. The vertical 



dashed line represents the finite jump in ks at the tangent bifurcation where the solution 7ij disappears. 



to have I-agents with a relatively low number of neighbours. As the infection is progressively suppressed 
by reconnection, the number of I-agents decreases and, accordingly, their contribution to the average 
number of neighbours per agent becomes less significant. At the bifurcation and beyond, S-agents must 
account for the whole average, so that ks returns to its value for f — 0, i.e. kg = k. 

The connectivity of I-agents, in turn, is a monotonically decreasing function of r, and reaches 
k\ = k/2X < k at the bifurcation. This implies that, even at the threshold of infection suppression, 
I-agents maintain a finite number of neighbours within the population. 

For A > 3/4, again, the connectivity kg associated with the solution rij initially increases with 

f, and attains a maximum. In the subsequent decay, however, it does not reach ks = k. In fact, nj 1 ^ 
disappears through a tangent bifurcation when it is still positive, so that the jump in the infection 
level is discontinuous. At the bifurcation, we find ks = 2k(2X + 3)/9 > k. The connectivity of I-agents 
decreases with r and, at the bifurcation, its value is independent of A: k\ = 2/3. 

From the viewpoint of the interplay of the epidemiological dynamics and the structure of the 
underlying network, the most interesting result of this analysis is the fact that the infection dies out 
even when infected agents keep a substantial connectivity with the rest of the population. In the cases 
illustrated in Fig. [4j for instance, infected agents preserve more than 60 % of their connections at the 
threshold where the infection level vanishes. In other words, as we had already verified for link deletion, 
reconnection needs not to completely isolate infected agents to suppress the infection. A moderate, 
partial isolation of the infected population is enough to asymptotically inhibit the endemic state. 



4 Heuristic description of infection suppression by link reconnection 

What mechanisms are at work when the infection is suppressed by reconnection, even when the con- 
nectivity of infected agents remains fairly high? To advance an answer to this question, it helps to 
consider a simpler dynamical system for the fraction of I-agents: 



n[ = -ni + [2A - f(l - rii) 2 ]ni(l - m). 



(20) 
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Fig. 5 Graphical solution of Eq. (|23|l . The dotted curve represents the right-hand side of the equation, and 
full curves are possible graphs of the left-hand side. Dots stand at their intersections. The arrows illustrate how 
the graphs may change upon the variation of parameters, in the cases of a transcritical bifurcation (A) and of 
a tangent bifurcation (B). 



The right-hand side of this equation is just a rearrangement of that of Eq. (fl4|) . The equilibria of Eq. 
([20ll are thus identical to the equilibria for m in Eqs. Moreover, their stability properties are also 
the same as in our original system. It is important to understand, however, that (|20|) and are not 
equivalent: they merely share the same equilibrium behaviour in what regards the fraction of I-agents. 

We immediately see that Eq. (|20|) can be put in the form of the standard mean-field equation 
for a SIS process if we introduce the effective infection probability 

A cff = 7 [2A - f(l - rn) 2 ] = (1 - r)kp - r(l - m) 2 . (21) 

In Eq. |(2J) the threshold of infection suppression, where the trivial equilibrium changes its stability, is 
given by A = 7. Imposing this same condition to A c ff, we find f = 2A — 1. But this is precisely the 
suppression threshold in the system with reconnections. Therefore, with respect to the stabilization 
of the trivial equilibrium, the system (jlip is effectively equivalent to the standard SIS model with 

infectivity A e fi- The transcritical bifurcation of Eqs. (flT]) . where = becomes stable, can be 
interpreted as a kind of continuation for r ^ of the transcritical bifurcation of the SIS model without 
reconnection events. 

The interpretation of the tangent bifurcation where the endemic state disappears at a positive 

value of Uj , for A > 3/4, is less direct. It can however be argued that the presence of such a bifurca- 
tion, together with the transcritical bifurcation which stabilizes the trivial equilibrium, constitutes the 
most generic critical behaviour expected for an epidemiological model like Eq. ^ when the infection 
probability depends on the density of I-agents: 

n[ = — 7?ii + A(nj)ni(l — nj). (22) 

Besides the trivial equilibrium, this equation has fixed points at the solutions of 

7- 1 A(nf) = (l-nJ)- 1 - (23) 

Figure [5] illustrates graphically two representative situations. The dotted curve is the graph of the 
right-hand side of Eq. (|2"3"1) as a function of n\. If the graph of 7~ 1 A(nj ) has a single intersection with 
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the dotted curve (A) and if, upon variation of parameters in the infection probability, the graph varies 
as indicated by the arrow, the intersection crosses n\ = and a transcritical bifurcation takes place. 
The standard SIS model, in which A is constant, is an example of this situation. More generally, the 
graph of 7 _1 A(nj ) may have two (B) or more intersections with the dotted curve. When the parameters 
change, it is still possible than one of the intersections becomes involved in a transcritical bifurcation 
crossing n\ — 0, as in situation A. Now, however, it may well be the case that two intersections 
approach each other, and eventually collapse and disappear, as in B. In this case, Eq. (f22|) undergoes 
a tangent bifurcation, as found to happen in our system 

In summary, the above discussion shows that the suppression of the endemic state as a result of 
reconnection events can be heuristically understood in terms of the critical behaviour of a standard 
SIS model with an effective infectivity, which depends on both the reconnection probability and on 
the fraction of infected agents. The transcritical bifurcation which stabilizes the state where the infec- 
tion is completely inhibited is interpreted as a continuation of a similar transition in the absence of 
reconnection. In turn, the tangent bifurcation -which, for large infectivities, suppresses the infection 
as the reconnection probability grows- is a generic phenomenon in SIS models with density-dependent 
infectivity [T5] . 



5 Conclusion 

In this paper, we have studied a model for an SIS epidemiological process in a population of agents 
on a network, where contagion can occur along the network links. The network coevolves with the 
population as the infection progresses: as a response to risk perception, susceptible agents can decide 
to break links with their infected peers. In the first version of the model, broken links are permanently 
removed from the network. For any positive probability of link removal, the infection is found to 
asymptotically die out. During the process, a fraction of network links is deleted, so that the social 
structure is degraded in the long-time limit. As expected, the fraction of remaining links decreases as 
the removal probability and the infectivity grow. 

In the second version, a susceptible agent who has broken a link with an infected agent reconnects it 
to a randomly chosen member of the remaining population. In this case, whether the infection persists 
or is asymptotically suppressed depends on the reconnection probability. Suppression of the endemic 
state does not require full isolation of the infected population. On the contrary, it can be achieved 
while each infected agent preserves a substantial part of the links with the rest of the population. 

Reconnection of network links introduces new dynamical features with respect to the standard SIS 
process. In particular, for sufficiently high reconnection probabilities, the continuous transition associ- 
ated with infection suppression is replaced by a discontinuous tangent bifurcation, where the infected 
fraction of the population drops abruptly as the relevant parameters are changed. The appearance 
of this new critical phenomenon is accompanied by the creation of a bistability regime, where the 
infection can either persist or die out depending on the initial fraction of infected agents. 

These results are remarkably robust under variations of the dynamical rules of the model. Here, for 
instance, we have assumed that link reconnection occurs, at each time step, before infection events, 
so that contagion can only take place from those infected agents who have retained their links. If 
this ordering is altered, the position of the bifurcation -and, consequently, the threshold of infection 
suppression- change, but the overall qualitative picture of Fig. [3J is not modified. The same holds 
if reconnected links are kept, with a certain probability, by infected agents, instead of always being 
susceptible agents which retain broken contacts. 

As shown in Section [31 the SIS model with link reconnections can be analytically solved in equilib- 
rium. To a large extent, the possibility of explicitly writing the stationary solution for the density of 
infected agents and the expression for critical lines in the parameter space is a direct consequence of 
the approximation done in Section[2l just before Eqs. (O, on the average number of infected neighbours 
of susceptible agents. As stated there, the approximation is valid as long as the number of infected 
agents remains high. On the other hand, we have pushed the solution to describe also the regime 
where the infection dies out. The question thus arises on whether our analytical solution still gives 
a reasonable description of the epidemiological dynamics. To advance an answer, we have performed 
agent-based numerical simulations of the SIS process on the evolving network. Results exhibit quanti- 
tative differences with the analytical solution, especially in the prediction of the critical points. Those 
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are precisely the zones where the analytical description is expected to fail, because at those points 
the infection is suppressed and the approximation breaks down. Qualitatively, however, the numerical 
and analytical results for the critical behaviour are the same. Simulations confirm the prolongation of 
the transcritical bifurcation for non-zero reconnection probabilities, and the appearance of the tangent 
bifurcation, along with the bistability regime, as reconnections become more frequent. 

Numerical simulations automatically incorporate further dynamical elements which are not present 
in a mean-field analytical description, such as the creation of correlations between the relative positions 
of susceptible or infected agents and the effects of heterogeneity in the distribution of network links. 
Extensive numerical results, as well as an improved analytical approach able to take into account 
correlations induced by the spatial structure of the network [15j . will be presented in a forthcoming 
paper [14]. 
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